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The  elecrrcnsrafic  Stability  properties  of  a rotating,  charge-neutralized 
P-layer  are  investigated  within  the  framework  of  a hybrid  (Vlasov-fluid) 
model  in  which  the  layer  ions  are  described  by  the  Vlasov  equation,  and  the 
layer  electrons  and  the  uniform  background  plasma  are  described  as  macroscopic, 
cold  fluids.  It  is  assumed  that  the  P-layer  is  thin,  with  radial  thickness 

Ac*  -v 

(2a)  much  smaller  than  the  mean  radius  (R«) , and  that  v <<  1,  where  v is 

d 

Budker's  parameter  for  the  layer  ions.  Electrostatic  stability  properties  are 
calculated  for  perturbations  about  a weakly  diamagnetic  P-layer  with  rectangular 
density  prof  ile* '■described  by  the  equilibrium  distribution  function 
f^-(ntRQ/2Trmi)6  [H-V^Pz-mi(VQ-V^)/2]<S(P0-Po) , where  H is  the  energy,  P0  is 
the  canonical  angular  momentum,,  Pz  is  the  axial  canonical  momentum,  and  n^, 

Rq,  Vz,  V0,  and  Pq  are  constants.  J The  stability  analysis  is  carried  out 
including  the  effects  of  a uniform  background  plasma,  and  weak  self  magnetic 

fields.  Although  a slow  rotational  P-layer  (P^>0)  is  found  to  be  stable, 

S- — j..t-  <j> 

it  is  shown  that  a fast  rotational  P-layer  (P^  < 0)  is  unstable  for  sufficiently 

/-"2 

high  background  plasma  density  '(u)p»w  )j.  T^e  typical  instability  growth 

rate  is  a substantial  fraction  of  the  ion  cyclotron  frequency.^ 
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I.  INTRODUCTION 

The  generation  ^ and  application^  ^ of  intense  ion  beams  has  been 
the  subject  of  several  recent  investigations.  One  application  of  con- 
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siderable  interest  is  the  formation  of  field-reversed  ion  layers  and  rings, 

which  can  provide  the  magnetic  confinement  geometry  for  fusion  plasmas. 

Such  layers  and  rings  are  likely  subject  to  various  macro-  and  micro- 

instabilities.'^  ^ For  example,  the  low  frequency  stability  properties 

( | w | <<wc^)  of  an  ion  layer  immersed  in  a background  plasma  have  recently 

been  investigated  within  the  framework  of  a kinetic  energy  principle  that 

incorporates  the  effects  of  large  ion  orbits. ^ In  the  present  analysis, 

allowing  for  perturbations  of  moderate  frequency  ( | wl^w^) , we  examine 
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the  equilibrium  and  negative-mass  stability  properties  of  a weakly 

diamagnetic,  charge-neutralized  proton  layer  (P-layer)  within  the 
framework  of  a hybrid  (Vlasov-fluid)  model.  The  layer  electrons  and 
background  plasma  electrons  and  ions  are  described  as  macroscopic, 
cold  fluids  immersed  in  a uniform  axial  magnetic  field  B^e^.  However,  to 
correctly  include  the  influence  of  layer  ion  dynamics  on  stability  behavior, 
we  adopt  a fully  kinetic  model  in  which  the  layer  ions  are  described 
by  the  Vlasov  equation. 

The  present  analysis  is  carried  out  for  an  infinitely  long  P-layer 
aligned  parallel  to  a uniform  magnetic  field  B^z  (Fig.  1).  The  P-layer 
is  immersed  in  a uniform,  cylindrical  background  plasma  and  is  charge  neu- 
tralized by  extra  electrons.  We  assume  that  the  layer  is  thin  [Eq.  (1)1, 
i.e.,  the  radial  thickness  (2a)  of  the  layer  is  small  in  comparison 
with  the  mean  radius  R^.  It  is  also  assumed  that  v<<1,  where  v is 
Budker's  parameter  for  the  layer  ions.  Although  the  equilibrium  self 


3 


magnetic  field  bqq  (r^e+Boz^r^z  is  weak  in  absolute  intensity  for  v<<1, 
the  self-field  gradients  can  be  sufficiently  large  to  have  an  important 
influence  on  particle  trajectories,  and  hence  on  stability  behavior. 
Equilibrium  and  stability  properties  are  calculated  for  the  specific 
choice  of  ion  layer  distribution  function  [Eq.  (10)], 


fb(H.Pz-V 


nbR0 

211111, 


O O 

“-V*  - T (VV 


4<Vpo>  • 


where  H is  the  energy,  Pa  is  the  canonical  angular  momentum,  P is  the 

o z 

axial  canonical  momentum,  and  n,  , R_,  V , V.,  and  Pn  are  constants. 

D U Z U U 

Equilibrium  properties  are  examined  in  Sec.  II.  One  of  the  important 

features  of  the  equilibrium  analysis  is  that  the  equilibrium  distribution 

function  in  Eq.  (10)  corresponds  to  a sharp-boundary  density  profile 

[Eq.  (29)].  It  is  also  found  that  the  radial  betatron  frequency  of  the 

2 22  1/2 

layer  ions  is  given  by  (oo^^^+oo^B  ) , where  w is  the  ion  cyclotron 

frequency,  is  the  ion  layer  plasma  frequency,  and  £ is  the  ratio 
of  the  mean  ion  layer  speed  to  the  speed  of  light. 

The  electrostatic  stability  properties  of  the  layer-plasma  configuration 
are  investigated  in  Secs.  Ill  and  IV,  including  the  important  influence  of: 

(a)  equilibrium  self  magnetic  fields,  (b)  an  outer  cylindrical  conductor 
(Fig.  1),  and  (c)  ion  layer  kinetic  effects.  The  analysis  is  carried  out 
within  the  framework  of  the  linearized  Vlasov-fluid  and  Poisson  equations, 
assuming  that  perturbed  quantities  are  independent  of  axial  coordinate 
(3/3z**0).  Moreover,  in  Sec.  IV,  stability  properties  are  investigated  for 
eigenfrequencies  near  multiples  of  the  mean  P-layer  rotational  frequency, 
i.e.,  |w-A(jOq  | where  w is  the  complex  eigenf requency , j l is  the 

azimuthal  harmonic  number  and  is  the  mean  rotational  frequency  of  the 
P-layer.  Although  a slow  rotational  equilibrium  (PQ>0)  is  stable  [see 
discussion  following  Eq.  (55)],  it  is  shown  in  Sec.  IV  that  a fast  rotational 


4 


equilibrium  (PQ< 0)  is  unstable  for  sufficiently  high  background  plasma 


2 2 

density  (u)p>>uc^) . The  physical  mechanism  for  instability  is  similar  to  that 

12-15 

for  the  negative-mass  instability,  including  the  dielectric 

effects  of  the  background  plasma.  Introducing  the  dimensionless  parameters 

[Eq.  (60)] 


,.2  n ^ “ci  _ 2.2.  2 

y-U  -1)  n — , q-u>bB  /u  , 


where  n^  and  n^  are  the  layer  density  and  plasma  density,  respectively. 


we  find  that  [Eqs.  (62)  and  (63)] 


y<l,  and  (l+q)y  -(3q/2+l)y+q>0  , 


are  necessary  and  sufficient  conditions  for  instability.  The  instability 
condition  in  Eq.  (62)  is  valid  only  when  the  parameter  n^/n^  is  sufficiently 
small  (n^/n  <<  1).  Moreover,  the  system  is  most  unstable  when  [Eq.  (65)] 


(VnP)=(1+q)/a  -1)  , 


which  corresponds  to  y=l.  Evidently,  the  mean  motion  of  the  ion  layer 
2 

(3  ),  the  background  plasma  (n^) , and  equilibrium  self  field  effects  (q) , 
all  have  an  important  influence  on  stability  behavior  (Sec.  IV). 

Numerical  investigations  of  the  stability  properties  are  carried  out 
in  Sec.  IV.  Several  points  are  noteworthy  in  this  regard.  First,  the  in- 
stability growth  rate  increases  when  the  self-field  strength  (as  measured 
by  q)  is  increased.  Moreover,  the  system  is  stabilized  as  q approaches  zero. 
Second,  the  maximum  growth  rate  can  be  a substantial  fraction  of  the  ion 
cyclotron  frequency.  Third,  the  number  of  unstable  modes  increase  rapidly 
as  q increases.  Fourth,  the  range  of  n^/n^  corresponding  to  instability 
is  rapidly  reduced  when  the  azimuthal  harmonic  number  l is  increased  above 
1=2.  [The  fundamental  mode  (*•=!)  is  found  to  be  stable.] 


I 1 . THEORETICAL  MODEL  AND  EQUILIBRIUM  PROPERTIES 


A.  Theoretical  Model 


As  illustrated  in  Fig.  1,  the  equilibrium  configuration  consists 
of  a nonrelativist ic  P-layer  that  is  infinite  in  axial  extent  and  aligned 
parallel  to  a unilorm  applied  magnetic  field  The  P-layer  is  immersed 


in  a uniform,  cylindrical  background  plasma  (with  outer  radius  Rq) , and  is 


charge  neutralized  by  extra  electrons  with  density  profile  identical 
to  the  layer  ions.  The  mean  radius  and  radial  thickness  of  the  P-layer  are 
denoted  by  R^  and  la,  respectively.  The  radius  of  the  cylindrical  conducting 


wall  is  denoted  by  R . The  mean  motion  of  the  P-layer  is  in  the  azimuthal 

c 


and  axial  directions,  and  the  applied  magnetic  field  provides  radial  confine- 
ment of  the  layer  ions.  For  simplicity,  we  assume  that  the  plasma  and  layer 
ions  are  singly  charged.  As  shown  in  Fig.  1,  cylindrical  polar  coordinates 
(r,t?,z)  are  introduced,  3nd  the  following  are  the  main  assumptions  pertaining 
to  the  equilibrium  configuration: 

(a)  Equilibrium  properties  are  independent  of  z (9/9 z=0)  and  azimuthally 
symmetric  (9/90=0)  about  the  z-axis. 

(b)  The  radial  thickness  of  the  P-layer  is  much  smaller  than  its  major 
radius,  i . e . , 


a « Rq  . 


(1) 


(e)  It  is  further  assumed  that 


2 2 

v =N,  e /m . c~  < < 1 
b l 


(2) 


2 2 

where  v =N^e  /m^c  is  Budker's  parameter, 


R 


"Ho 


NY  — 2tt I C dr  r n°(r) 

D 


(3) 


is  the  number  of  layer  ions  per  unit  axial  length  of  the  P-layer,  nU(r) 

b 


1 


1 


6 


is  the  layer  ion  density,  c is  the  speed  of  light  in  vacuo,  and  e and  iik 
are  charge  and  mass,  respectively,  of  a layer  ion.  The  inequality  in  Eq . 
(2)  implies  that  the  intensity  of  the  self  magnetic  field  is  much  less  than 
the  applied  field  BQ  (i.e.,  AB0/BQ«1).  However,  the  self-field  gradients 
over  the  narrow  radial  dimension  of  the  layer  can  be  sufficiently  strong 
to  have  a large  influence  on  the  layer  ion  trajectories,  and  hence  on 
stability  behavior  (Sec.  III. A). 

In  the  present  analysis,  the  layer  electrons  and 
background  plasma  are  treated  as  cold  (T  -*•())  fluids  immersed  in  a uniform 
axial  magnetic  field  B^^.  Here,  j=e',  e,  and  i represent  the  layer 
electrons,  plasma  electrons,  and  plasma  ions,  respectively.  Within  the 
context  of  the  electrostatic  approximation  (V*£=0) . the  equation 
of  motion  and  the  continuity  equation  for  each  fluid  component  (j=e',  e,i) 
can  be  expressed  as 


(—  + v • v)v 
V 3t  vfy 


(4) 


+v* 

'V 


<W° 


(5) 


where  J£(;jc,  t ) =— V<(>  (^c,  t ) is  the  electric  field,  Uj(^c,t)  is  the  density,  ^(}c,t) 
is  the  mean  velocity,  and  and  rrn  are  the  charge  and  mass,  respectivel; 

of  a particle  of  species  j.  In  Eq.  (4),  the  electrostatic  analysis  is 
consistent  to  the  low  beta  approximation. 

To  include  the  influence  of  layer  ion  dynamics  on  stability  behavior, 
we  adopt  a fully  kinetic  model  in  which  the  ion  layer  distribution 
function  f,  (x,v,t)  evolves  according  to  the  Vlasov  equation 

D v 'w 


(3  . S . « / „ . rW\  a ) 

« + k • » + ;rr'* + — — j ■ « 

*\j  1 'XjJ 


fb(^»t)=0  • 


(6) 
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In  Eqs.  (4)-(6),  the  electrostatic  potential  <t>(£,t)  is  determined  self- 
cons  Lstently  from  Poisson's  equation 


V di=-4ir 

O . T 


1JP  fb  + l en($,t)  . 

j=e',e,iJ  J > 


(7) 


Equations  (4 ) — ( 7 ) constitute  a closed  description  of  the  system  and  form 
the  theoretical  basis  for  the  subsequent  analysis. 


B . Equilibrium  Properties 

For  azimuthally  symmetric  equilibrium  profiles  characterized  by 


nj(r)  and  V^(x)=v9g  , j=e',e,i,  it  is  straightforward  to  show  from  Eq.  (5) 


that  the  functional  form  of  the  electron  layer  and  background  plasma  density 


0, 


profiles  n^ (r)  can  be  specified  arbitrarily.  Moreover,  from  equilibrium 


charge  neutrality,  the  equilibrium  radial  electric  field  vanishes,  i.e., 


E^(r)=-3<j9(r)/9r=0.  It  follows  from  Eq.  (4)  that  equilibrium  force  balance 


in  the  radial  direction  can  be  expressed  as 


U).(w.+e,oo  ,)=0  , 

1 J j cj 


(8) 


where  e =sgne  , and  io  (r)=V;  (r)/r  and  w =eBQ/m.c  are  the  angular  velocity 
J 3 J J 7 


and  cyclotron  frequency,  respectively.  In  obtaining  Eq.  (8),  we  have 
approximated  ^0(^)=BQ^z  consistent  with  Eq . (2).  In  general  we  note 


from  Eq.  (8)  that  there  are  two  allowed  equilibrium  values  of  Wj(r).  Throughout 


the  subsequent  analysis  we  assume  that  the  layer  electrons  and  background 
plasma  components  are  rotating  in  the  slow  rotational  mode  with 


Wj(r)=0  , j=e\ e, i . 


(9) 


For  the  layer  ions,  any  distribution  function  f,  (x,v)  that  is  a function 

b ^ 


only  of  the  single-particle  constants  of  the  motion  in  the  equilibrium 
fields  is  a solution  to  the  steady-state  (9/9t=0)  ion  Vlasov  equation. 
For  present  purpose,  we  assume  an  ion  layer  equilibrium  described  by 


J 


8 


fb°(H,Pz,P  ) 
b z 0 27Tm, 


i 2 2 

Vzpz  - T < VV 


6(pe_po)’ 


(10) 


where  V z and  VQ  are  constants,  nb  is  the  layer  density  at  the  equilibrium 
radius  r=RQ,  H is  the  total  energy, 


„ _ 1 ,2.2.2. 

H " 2^7  (pr+p6+pz)  ’ 
1 


P is  the  canonical  angular  momentum, 
0 


pe  " r 


m . 

pe  + Tuc±r  +tAeS(r> 


and  P is  the  axial  canonical  momentum, 
z ’ 


(ID 


(12) 


P 

z 


+ 


- A3 * S 
c z 


(r)  . 


(13) 


s s 

Here,  Ay(r)  and  Az(r)  are  the  6-  and  z-components  of  vector  potential  for 

the  axial  and  azimuthal  self  magnetic  fields.  Without  loss  of  generality, 

we  assume  that  the  vector  potentials  A„(r)  and  As(r)  vanish  at  r=R„. 

o z u 

In  Eqs.  (11)-(13),  lower  case  jj>=nuv  denotes  mechanical  momentum. 

For  the  choice  of  ion  distribution  function  in  Eq.  (10),  the  equilibrium 
vector  potentials  are  to  be  calculated  self-consistently  from  the  steady- 

g 

state  Maxwell  equations.  The  0-  and  z-components  of  the  Maxwe11 

equation  can  be  expressed  as 


3 13  s,  , 

"5 t — rA„  (r)  = 

3r  r 3r  0 


4ire  0.  . 0,  . 

— Vr)Vr) 


“r  fr  r ST  Az(r)  = 


4ire  ,3  ,0,„  „ „ \ 

— a P ve  fb<H.Pz,Pe)  , 


Aire  0,  .,.0,  . 
— nb(r)Vz(r) 


(14) 


d3p  vz  fb(H,Pz,pe)  • 


where  the  local  ion  layer  density  n^(r)  is  defined  by 
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n£(r)=Jd3p  fJ(H,Pz,Pe)  . (15) 

Here,  v =p  /m  , v =p  /m, , and  V^(r)  and  V^(r)  are  the  mean  azimuthal 
WW1ZZ1  u z 

and  axial  velocities  of  an  ion  layer  fluid  element. 

f 3 f°° 

Substituting  Eq . (10)  into  Eq . (15)  and  representing  d p=2tt  dp  dp^p^, 

2 2 2 u 

where  p^-p^+Cp^-nK  V^)  > we  find  that  the  ion  layer  density  profile  is  given 


Of  X 
nb(r)  = 


r<R1  , 

Rx<r<R2  , 

R <r<R  , 
l c 


where  R^  and  R^  are  the  extremes  of  the  interval  on  which  the  inequality 


(r)>0 


is  satisfied.  [That  is,  R^  and  R2  are  determined  from  ^(R^)=\p(R^)=0. ] 
In  Eq.  (17),  the  envelope  function  <j/(r)  is  defined  by 


m . „ , fPn  m z 

. l „2  1 0 i e .s,  , . e „ .S/  , 

) = -r~  V - — — — id  .r A.  (r)  + — V A (r)  . 

20  2m.r  2cx  c0  czz 

l ' 


It  is  evident  from  Eqs.  (10)  and  (18)  that 


m , « « 

H - V.  - f (V0-Vz> 


= [p^+(p  -m,V  )2H(r) 

2nu  r z i z 


for  Pfi=Pg.  That  is,  t|i(r)  is  the  (r-z)  kinetic  energy  of  a layer  ion 
in  a frame  of  reference  moving  with  axial  velocity  V . 


Thus  far,  RQ  has  been  introduced  in  the  analysis  as  an  unspecified 
constant  parameter  in  Eq.  (10).  Without  loss  of  generality,  we  now  choose 


10 


Rq  to  correspond  to  that  radius  where  ip(r)  passes  through  a maximum  in  the 
interval  R1<r<R?,  i.e., 


fe  *(r) 

l3r 


= 0 . 


(19) 


r=Rn 


Substituting  Eq.  (18)  into  Eq.  (19),  and  making  use  of  £q(;t)=Bq6^(+B® 
-■O/Br)  A^(r)^g+(l/r)  (3/3r)  [rAp  (r)  we  obtain 


we+'uCi+[eB0z (Ro) /mic  1 }tVVzeBoe (Ro) /micRo=0  ’ 


(20) 


where  a)0=Vg (R^) /Rq  is  the  angular  velocity  of  an  ion  layer  fluid  element  at 

+ 

t=Rq,  Solving  Eq.  (20)  for  u>g  gives  Wq=u)0,  where 


w„--w  . , 

0 Cl 


(21) 


corresponding  to  a fast  rotational  equilibrium  with 


1 2 

Pn  = - o’  m.Rnw  .<0 
0 2 i 0 ci 


(22) 


Similarly, 


VWV'Vo- 


(23) 


corresponding  to  a slow  rotational  equilibrium  with 


1 2 

Pn  = m.Rnw  .>0  . 
0 2 i 0 ci 


(24) 


In  obtaining  Eqs . (21)  and  (23),  use  has  been  made  of  Eq . (2),  which  implies 

weak  self  fields  with  | B® I , IbS|<<B.,.  Note  that  there  are  two  classes  of 

1 0 1 1 z 1 0 

equilibria  [Eqs.  (21)-(24)],  namely,  a fast  rotational  equilibrium  with 
Pq<0,  and  a slow  rotational  equilibrium  with  Pq>0. 

A closed  analytic  determination  of  R^  and  R£  from  the  zeros  of  Eq.  (18) 
is  not  generally  tractable  except  for  a thin  layer  [Eq.  (1)].  We  now  spec- 


ialize  to  the  case  of  a thin  P-layer  with  (R  -R^/R^d  and  Taylor  expand 
Eq.  (18)  about  r=RQ>  i.e., 

'f(r)  = \ mitVo“(Po/miR0_u>ciR0/2)21  ~1  mi“r(r-Ro)2+--- 


where 


2 _ 1 /3  > 

“r  " m.  ( 2 

x r 


2 2e  2 
=0)  ,+to,  8 

cl  b 


1 3—1  S 

is  the  radial  betatron  frequency-squared  of  a layer  ion.  In  Eq . (26) 

2.2.  ? 

“b-^6  nb'mic  ls  the  i°n  plasma  frequency-squared,  and  8 is  defined  by 


uu  c is  tn 
2_„02/n  ,,22,  2 


^ \J  £ / / / TOO 

8 =V0  (R0)/c  +vz/c  '(p0/miRoc_UciR0/2c)  +vz/c  • Xt  is  important  to 


recognize  that  the  term  in  Eq.  (26)  is  directly  related  to  the  self- 

field gradients  at  ^=Rq.  In  particular,  making  use  of  Eq.  (14),  it  is 

2 2 

straightforward  to  show  that  can  he  exPressed  as 


V - 

co,  6 = 

D 


2 c2  ^9B0z\  /I  3_  s \2 

“ci  22  (gr  ) (r  Jr  rB0e ) ' 

wKBn  IV  < r=R- 


Defining  the  half  thickness  of  the  P-layer  by 


a=tV0-(P0/miVU)ciR0/2)211/2/“I 


and  substituting  Eq.  (25)  into  Eq . (17),  we  readily  find  R^=RQ-a  and  R^= 
RQ+a,  where  a is  defined  in  Eq.  (27).  Finally,  we  evaluate  the  angular 

velocity  u>~  [Eq.  (23)]  for  a slow  rotational  equilibrium.  Making  use  of 

9 

Eq.  (14)  to  eliminate  Bq0(Ro),  we  obtain  from  Eq.  (23) 


we=Vwci(Vz/Vci)  ’ 


2 2 

where  v=Nbe  /true  is  Budker’s  parameter  defined  in  Eq.  (2),  and  Nb=4nnbR0a 
is  the  number  of  layer  ions  per  unit  axial  length. 
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For  a thin  layer,  the  ion  density  profile  in  Eq . (16)  can  be  approximated 
by 

(V  lr"Rol<a  • 

n°(r)=  J (29) 

( 0 , | t-Rq | >a  . 

Assuming  a space-charge  neutralized  layer,  the  layer  electron  density  profile 
ne'(r)  is  identical  to  Eq . (29).  Moreover,  the  background  plasma  has  a 
uniform  density  profile  with 

nj(r)=V  0<r<Rc>  (30) 

where  n^  is  constant.  We  now  investigate  electrostatic  stability  properties 
for  perturbations  about  the  plasma-layer  equilibrium  characterized  by 
Eqs.  (9),  (10),  (25),  (29),  and  (30). 


j 
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III.  ELECTROSTATIC  STABILITY  PROPERTIES 


A.  Eigenvalue  Equation 

In  this  section,  we  linearize  Eqs.  (4 ) — ( 7)  assuming  electrostatic 
perturbation  about  the  equilibrium  described  by  Eqs.  (9),  (10),  (25),  (29), 
and  (30).  The  present  analysis  assumes  flute  perturbations  with  3/3z=0, 
so  that  all  perturbations  have  spatial  dependence  only  on  the  perpendicular 
variable  x1=(r,0).  In  the  electrostatic  approximation,  the  perturbed 
electric  field  is  6E(x,t)=-V15<t>  (^,  t) , and  Eqs.  (4)-(7)  can  be  linearized 
to  give 

no 

Ft  Sni  + 7 h <rnj5vlt)  + r1  le  5v3e'°  ' 


s'V-'iW-ilfc1'  • 


(31) 


ft  svje  + hW  * s*  ■ 

for  the  layer  electrons  and  plasma  electrons  and  ions  (j=e\  e,i),  and 

3 . L0, ...  . * *W|  . 3_1  . , . 

3^/8fb(«A«*t) 


{fr 


+ K 


9£± 


+ e 


0 , v ^ ^ 

« + -7 


(32) 


-e[V<K&,t)]  * fb(^i’«) 

for  the  layer  ions  (j=b).  In  obtaining  Eq . (31),  we  have  approximated 
cons^stent  wit'1  Eq.  (2).  The  linearized  Poisson  equation  is 

(f  h r Tr + 7 p>) s*  - • 

where  use  has  been  made  of  Eq . (9).  In  Eqs.  (31)-(33),  ^j(^j.*t)  and 
6nj(^A,t)  are  the  perturbed  fluid  velocity  and  density,  6f^(^c1,j^,t)  is 
the  perturbed  ion  layer  distribution  function,  and  ^°(^)  =-(3i)>0/3r)^  =0 
is  the  equilibrium  radial  electric  field. 


(33) 
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To  simplify  the  right-hand  side  of  Eq.  (32),  use  is  made  of  3P Q/3£= 
reQ  and  3U/3j>=v-Vz£z,  where  l^H-V^P^.  Here  and  jg  are  unit  vectors  in 
the  9 and  z directions,  respectively.  We  express  perturbed  quantities  as 


6ij;(^±,t)=ii)Jl(r)exp{i(«.b-wt)  },  lmw>0, 


(34) 


in  Eq.  (32),  and  integrate  from  t'=-°°  to  t'=t,  using  the  method  of 
characteristics.  In  Eq.  (34),  u is  the  complex  eigenf requency , and  £ is 
the  azimuthal  harmonic  number.  Neglecting  initial  perturbations, 
and  noting  that  3f [V 3U  and  3f^/3P  are  constant  (independent  of  t')  along 
particle  trajectories  in  the  equilibrium  field  configuration,  the  perturbed 
ion  distribution  function  can  be  expressed  as 


3f 


= e nr  ^(r)  + ie r 


5 

3U 


+ l 


9Pn 


(35) 


where  the  orbit  integral  I is  defined  by 


(r ' ) exp  [-iwT+iJi  (6  ' -0  ) ] . 


(36) 


Here  T=t'-t,  and  the  particle  trajectories  ^'(t')  and  v'(t')  satisfy 
d^'/dt'=v'  and  m^dv'  /dt'=ev'x^(^c')/c,  with  "initial"  conditions  ^'(t^t)*8^ 
and  (t'=t )=£. 

The  evaluation  of  the  orbit  integral  in  Eq.  (36)  is  generally  complicated. 
However,  for  present  purposes,  we  consider  low-frequency  perturbations 
satisfying 


2 2 

Iw-S-ioJ  <<w  , 

1 9'  r 

SLa/Rr.<<w  /w  . , 
0 r ci 


(37) 


where  a is  the  half-thickness  of  the  layer  defined  in  Eq.  (27).  Within 
the  context  of  Eq.  (37).  it  is  valid  to  approximate15 
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, , .w  ,»  W 6P0  W0/3^\ 

Vr  )"*t(R0)  2 \ 3r  / _ 


U u r-R_ 

r 0 


in  Eq.  (36) , 


e"6  + l“»'^/T 

In  Eq.  (38),  6P  =P  -P_ , and 

t)  o U 


U()  2P0/mR()-+Ucl  , 


2,2.  22.2 

y-a)c/wr-i=-WjjB  /ur  < 


where  is  the  radial  betatron  frequency  defined  in  Eq.  (26),  (+) 

refers  to  fast  and  slow  rotational  equilibria,  and  a>  is  the  angular 

0 

velocity  of  an  ion  Layer  fluid  element  at  r=R„ [u>  =uj  =-w  for  a fast 

0 0 6 ci 

rotational  equilibrium,  and  w =w  = vw  .(V  /R.w  .)2  for  a slow  rotational 

6 8 ci  z 0 ci 

equilibrium,  o.f.,  Eqs.  (21)  and  (28)].  Substituting  Eq.  (38)  into 

Eq.  (36),  and  approximating  (v0/r)p  _p  =a)^-w0(r-RQ) /Rq,  Eq.  (35)  can 

0 0 

be  expressed  as 


Vr = e 


3fb  f*  - *“0(r“V  . 

3U  <r)~^Jl<:R0)  “ (a)-Lai  )R  W 

0 u 


ef^  fw*.  ( 3<t>." 


mR  (uj-Huj  ) R,,(ao  Loo.)  2 V 9r 

•I  6 U 0 U)  ' 


Wo  further  assume  that  the  ion  layer  density  satisfies 


n <n 
b°-  p 


which  will  assure  the  validity  of  Eq.  (37)  in  the  subsequent  analysis. 
To  simplify  notation,  we  also  introduce  the  dielectric  function 
of  the  background  plasma 
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01 


e(oj) 


1 - 


= 1 + 


JH. 


0) 


2 2 2 2 
0)  -U)  0)  -0)  . 

ce  ci 


(42) 


0J 


0) 


o)  oi  . 2 2 

ce  ci  to  -oi  . 

ci 


2 2 

where  use  has  been  made  of  Eq.  (37),  and  (ji^=4ire  n^/iiK  is  the  background 
ion  plasma  frequency-squared.  Substituting  Eqs.  (34)  and  (40)  into 
Eqs . (31)-(33),  and  making  use  of  Eqs.  (37)  and  (42),  it  is  straight- 
forward to  express  the  eigenvalue  equation  as 


1 3_ 

r 3r 


re(o))—  ^(r) 


2 ^ / 03 

K E(u>)i(r)  = [6(r-R,)-6(r-R1)l 

l r wo)  , 2 1 

r 


R.0) 


0 


ci 
r-R 


+ iU(r)-^(R0)]  - ( } 


o : 


a 2 
0) 

r 


(R, 


0' 


0 


41 


6(r-R1) 


6 

+ — 


(r-R2)  ) 2 R0 

r 1 b r 


>[(R2-r)(r-R1)] 


dii2kR0) 


(43) 


2 2 
Ro(o)-Au)0) 


Ho) , 


R„o) 

0 r 


1 

{ 3r  jR 

(w-Hu)  ) 
n 0 

2 2 * A 

where  w =4Tre  n, /m . , <f>  (r)  H<j>  (r)  , and  ®(x)  is  the  Heaviside  step  function 
b b i x. 

defined  by 


(x)  = 


1 , 

x>0 

0 , 

x<0 

(44) 


In  obtaining  Eq . (43),  we  have  neglected  the  additional  electron  layer  contribu- 
2 2 2 

tion,  -oi  i(r)/[oj  -o)  . (r)],  to  the  dielectric  function  e(o>)  on  the  left-hand 

pe  ce 

side  of  Eq.  (43).  For  the  low  frequencies  considered  here,  the  corresponding 


term  is  of  order  (n,/n  ) (m  /m.)  smailer  than  the  plasma  ion  contribution  to  e(oi). 
bp  e i 

It  should  be  noted  that  the  terms  in  Eq.  (43)  proportional  to  6(r-R^) 

and  {(r-R^)  correspond  to  surface-charge  perturbations  on  the  inner  and 

outer  boundaries  of  the  P-layer.  These  terms  have  a decisive  influence  on 

2 2 

stability  behavior  when  w.  /w  >>1.  We  further  note  that  the  final  term 

b ci 


17 


on  the  right-hand  side  of  Eq.  (43)  is  proportional  to  «[  O^-r) (r-Rj) ] and 
corresponds  to  a body-charge  perturbation.  The  general  form  of  Eq . (43) 
is  similar  to  the  eigenvalue  equation  obtained  by  Davidson  et  al.^  for  a 
nonrelativistic  E-layer,  and  similar  techniques  can  be  used  to  determine 
the  complex  eigenf requency  u>. 


B.  Dispersion  Relation 

We  now  investigate  the  stability  properties  predicted  by  Eq.  (43). 
Since  the  surface  terms  in  Eq.  (43)  are  nonzero  only  at  the  boundaries 
of  the  P-layer,  the  perturbed  potential  at  all  other  radial  points  satisfies 

f 0 , 0<r<R1  , 


(pi2$ (Rq)  *u0Rq  (^/8r)R0] 


EV 


(id-JlWg)' 


(w-£w0) 


Vr<R2  » 


l 0 , R2<r<Rc  » 


(45) 


where  R^=Rg-a  and  R2=RQ+a  are  the  boundaries  of  the  P-layer.  For  £>2, 
the  physically  acceptable  solution  to  Eq.  (45)  is 
f Ar^  , 0<r<Rj  , 


^ (r)  = < 


u “*■ 

Br  +cr 


2, 

“b/e  r 
t2-l  R0 


2- 

hi  4>(R0) 

[(.ii-iw  )2 


ici)0R0O<()/3r)R 


U)  (u>-£(D„) 

r 6 


R1<r<R2  , 


V -I 

^Dr  +Er  * 


R„<r<R 
2 c 
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where  A',  B',  C',  D',  and  E1  are  constants.  Since  Eq.  (43)  Is 
similar  in  general  form  to  Eq.  (88)  of  Ref.  15,  we  briefly  outline 
the  derivation  of  the  dispersion  relation.  [The  detailed  procedure 
is  similar  to  Ref.  15- ] 

We  introduce  the  abbreviated  notation 

Rn  “k 

Sl^  -SF-2  ’ 

ur 

2 

w,  !lwn 

s,U)  = - s , 

2 2 ((jj-e.a)„) 

we  9 


S,(w)  = — 1 - 
3 eww 


^2 

N (w)  = \i  — ~ > 

i -1  (w-f.w.) 

U 

2 

, (w,/e)iwn 

N.(w)  = - -J-  ^ . 

£ -1  w (w-f.w  ) 

r 9 


Three  constraints  relating  the  constants  A,  B,  C,  D,  and  E (or  A', 

B’ , C' , D' , and  E')  are  obtained  from  the  boundary  condition 

i(r=Rc)=0,  and  by  enforcing  continuity  of  4>(r)  at  r=R^  and  r=R2  • 

The  remaining  two  constraints  are  obtained  by  integrating  Eq.  (43) 

across  the  layer  boundaries  at  r=R^  and  at  r=R2»  After  some  tedious 

algebra,  we  obtain  the  dispersion  relation 

2 1 2 1 
|4i2gf+2US1+S3)+2il(S1-S3)  j + (2Hgf+S1+S3)(S1-S3)  1 ~ j 

/Ri\*  /R0\ 

x(l+N1+N2)  + (2igf+S1+S3)(x1-i3)  ^-R-  J - (2«.+S1-S3)(x2+X4)  ^ J 

n 

+ (S1-S3)  (x2-X4)  ~ [2*  (gf-l)  + (S1+S3)  (Xx+X3)  ] ^ r2  ^ 
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-R,\  C/R 


+2t,<2  (ic)  (id)  (1-1-S1+S3)<S1-S2)+(r:)  <l-W2»Sf«1+S3)(S1+S2)  -0, 


'2'  l “0 


for  i> 2.  In  Eq.  (47),  Xj (w)  is  defined  by 


Xl(w)  = r“  <1-i-s1+s3)N1-(S1+S2)(l+N2)  , 


X2(w)  = — (l-JH-2Jlgf+S1+S3)N1+(S1-S2)(l+N2) 


X3(w)=  jj-  (I-H-S1+S3)+S1+S2  2.N2  , 


X4(u>)=  — (l-fc+2«,gf+S1+S3)-S1+S2  £N2 


and  the  geometric  factor  g^  is  defined  by 


g-[l-(R,/R  )2SL]  1 . 
t l c 


Equation  (47)  is  the  desired  dispersion  relation  for  i>2.  Keep  in  mind 


that  £<<Rg/a  has  been  assumed  in  deriving  Eq.  (47)  [see  Eq.  (37)]. 


Substituting  R^R^-a  and  R^R^+a,  and  Taylor  expanding  the  left-hand 


side  of  Eq.  (47),  for  a<<RQ,  we  find  that  the  dispersion  relation  can 


be  approximated  by 


-2  f *+r3  -T-  £ ’ 0 • 

eo)  0 eco  . 0 

r ci 


where 


x=(o)-£a). ) /a)  . 

n Cl 


is  the  normalized  Doppler-shifted  eigenfrequency , I\  is  defined  by 


j 


2 
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a_ 

R, 


r^gf.w)  gf  (!  + 2 ) + 2 ( 2 . 

Eoj  / \ ew  / 0 

r / \ r ' 


r2(gf,io)=  2 + — 

ew 


fi.  (g_— 1)  + 

°f  EWU) 


Cl 


2eu 


2 * 


(52) 


r,(u))-n|i  + T 


2 

4 \ 
“b  \ 

CM  U| 

3 | 

/ 2 4 
f“b  % 

2 + 

2 4 J 

2 I 

2 2 

4 

eo) 

E 0)  / 

(a) 

\eo)  e 

03 

r 

r/ 

r 

\ r 

r 

and  (+)  refers  to  fast  and  slow  rotational  equilibria.  In  obtaining 
Eq.  (50),  use  has  been  made  of  the  definitions  in  Eqs.  (39),  (46),  and 
(48). 


The  preceding  analysis  pertains  to  i>2.  A similar  procedure  can  be 
followed  to  obtain  the  dispersion  relation  for  5.=1.  After  some  algebraic 
manipulation,  it  is  straightforward  to  show  in  the  thin-layer  approximation 
that  the  dispersion  relation  in  Eq.  (50)  is  also  valid  for  i=l.  In  this 
context,  Eq.  (50)  is  the  desired  dispersion  relation,  which  determines 
the  electrostatic  stability  properties  for  l<£«Rg/a. 
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IV.  STABILITY  ANALYSIS 


We  now  investigate  the  stability  properties  predicted  by  Eqs . (50)- 
(52).  The  stability  analysis  for  a slow  rotational  equilibrium 
proceeds  in  the  following  manner.  Since  the  angular  velocity  of  an  ion 
layer  fluid  element  for  a slow  rotational  equilibrium  is  much  less 
than  the  ion  cyclotron  frequency,  i.e.. 


U)  =vw  .(V  /R„w  .)  <<(D  . , 
0 ci  z 0 ci  ci 


it  is  readily  shown  from  Eq.  (42)  that  the  background  plasma  dielectric 
function  can  be  approximated  by 


e(i.w.)-l+(l+w  . / a)  )(w2/w2.)>0  . 
0 ci  ce  p ci 


Substituting  Eq.  (53)  into  Eq.  (52),  yields  the  inequalities 


rx(gf ,Jlu)e)>0,  r3(«.Wg)<o  . 


Within  the  context  of  Eqs.  (50),  (53),  and  (54),  it  is  straightforward 
to  show  that  the  system  is  stable  for  a slow  rotational  equilibrium. 

The  stability  analysis  for  a fast  rotational  P-layer  is  generally 
more  complicated.  We  therefore  restrict  the  analysis  to 

2 2 

<*>  »u  . . (55) 

p ci 

Making  use  of  n^n^  in  Eq-  (41),  it  is  straightforward  to  show  that 

lavr^/ar.r-eu2. |<i. 

b 2 13  ci1^ 

In  this  context,  the  dispersion  relation  can  be  approximated  by 


r x2+r  M 

13  2 R 

ew  . 0 

ci 


= 0 , 
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where  the  term  proportional  to  x has  been  neglected.  Moreover, 
making  use  of  the  inequality  in  Eq.  (37)  the  dielectric  function  in 
Eq.  (42)  can  be  approximated  by 


h 


(x2-2x)-1 


A-l  , 


: (w)  ~- 


(57) 


“ci/  f [x2-2Hx+(Jl2-l)  ]_1,  i> 2 , 


where  x=  (uj-£id„  ) /u)  . = (u)+i.u>  . )/w  .. 

0 ci  ci  ci 


A general  numerical  analysis  of  Eqs.  (56)  and  (57)  is  summarized 
later  in  this  section.  However,  for  eigenf requencies  satisfying 


w+Jlw  . <<  u . , 

Cl  Cl 


(58) 


Eq.  (57)  tan  be  further  simplified  to  give 

2 


e(u) - 


\ 

*) 


Cl 


~(2x)  1 , 

l-l  , 

(59) 

U2-l)-1  , 

l>2  , 

and  the  solution  to  Eq.  (56)  is  analytically  tractable.  Defining 


y=(^2-l) (nb/np) (uci/wr)2  , 


(60) 


2 2 2 2 2 2 
q=((Vu)ci)/a)ci=lJb6  /uci  » 


Eq.  (56)  can  be  expressed  as 

2U 


. 2 la 
(l-y)x  + — 

R0 


1 -(^)  y (q+l)y2-(|  q+ijy+q  =o  , (61) 


for  ?>2.  Here,  q is  a measure  of  the  strength  of  the  self-field 


2„2 


gradients  at  r=Rg.  [See  expression  for  u>  B following  Eq.  (26).] 


In  obtaining  Eq.  (61),  use  has  been  made  of  Eqs.  (1),  (26),  (39), 
(49),  and  (59).  Note  from  Eq . (60)  that  y is  positive  for  l>2. 


A parallel  analysis  for  1=1  results  in  a dispersion  relation  with 
stable  solutions  in  the  parameter  regime  of  physical  interest.  We 
therefore  restrict  discussion  to  mode  numbers  l>2. 

We  note  from  Eq.  (61)  that  (q+l)y2-(3q/2+l)y+q>q/2  for  y>l. 

Therefore,  the  necessary  and  sufficient  conditions  for  instability  can 
be  expressed  as 

y<l  , (62) 

and 

(l+q)y2-(|  q+l)y+q>0  , (63) 

for  l>2  and  | | <<  w It  is  straightforward  to  show  that  Eq. 

(63)  is  automatically  satisfied  for 

22  2 3/2 

q=wb6  /w^>2(2 -1)  /7=0. 52 . (64) 

For  the  case  q>0.52,  the  necessary  and  sufficient  condition  for 
instability  is  given  by  y<L  in  Eq.  (62).  It  should  be  noted  from  Eqs. 

(62)  and  (63)  that  the  stability  criteria  are  independent  of  the  location 
of  the  conducting  wall.  Moreover,  the  system  is  stable  when  the  self 
field  is  negligibly  small,  i.e.,  q<<l.  Finally,  we  note  that  the  maximum 
growth  rate  occurs  at  the  resonant  value 

(nb/np)  = (l+q)/(Jl2-l)  , (65) 

corresponding  to  y=l. 

Removing  the  restriction  in  Eq.  (58),  the  dispersion  relation  in 

Eq.  (56)  has  been  solved  numerically  by  substituting  Eq.  (57)  into 

Eq.  (56).  Figure  2 shows  a plot  of  the  normalized  growth  rate  Imu  /to  .=m./w  . 

c x i ci 

versus  n,/n  for  1= 2,  a/R  =0.05,  R„/R  =0.5,  and  several  values 
bp  0 0 c 
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of  q.  Several  points  are  noteworthy  in  Fig.  2.  First,  the  instability 

growth  rate  is  considerably  increased  when  the  strength  of  the  self 

magnetic  field  (as  measured  by  q)  is  increased.  For  example,  the 

maximum  growth  rate  (w.)  =0.114  w . for  q=0.5  occurs  at  n,  / n =0.63, 

l m ci  bp 

while  the  maximum  growth  rate  (w.)  =0.25  w . for  q=2  occurs  at  n,  /n  = 

l m ci  bp 

0.65.  We  also  note  that  the  value  of  n^/n^  corresponding  to  maximum 

growth  increases  as  the  self  magnetic  field  is  increased  [see  also 

Eq.  (65)].  Second,  the  maximum  growth  rate  can  be  a substantial 

fraction  of  the  ion  cyclotron  frequency  Third,  even  for  weak 

self-field  effects  (q<l),  the  range  of  n^/n  corresponding  to  instability 

is  considerably  extended  beyond  the  range  of  n^/n^  predicted  by  Eq.  (62). 

For  example,  for  q=0.5,  Eq.  (62)  predicts  that  instability  occurs  only 

for  n^/np<0.5.  On  the  other  hand,  from  Fig.  2,  we  note  that  the 

system  is  unstable  for  values  of  n,  /n  up  to  n./n  =1.6.  However,  the 

b p b p 

absolute  maximum  growth  rate  does  occur  very  near  to  the  resonant  value 

of  n^/np  (for  example,  n^/n^O-S,  for  q=0.5)  [Eq.  (65)].  Fourth,  the  growth 

rate  curve  has  more  than  one  maximum  for  q<0.52,  which  is  predicted  by 

Eq.  (63).  [See  the  curves  corresponding  to  q=0.25  and  q=0.5  in  Fig.  2.] 

We  also  emphasize  that  the  range  of  n^/n^  corresponding  to  instability 

is  rapidly  reduced  when  the  azimuthal  harmonic  number  £ is  increased. 

Of  considerable  interest  for  experimental  application  is  t.ie 

stability  behavior  for  specified  values  of  n^/n^,  q.  and  R^/R^. 

Typical  results  are  shown  in  Fig.  3 where  Imw/w^  is  plotted  versus 

mode  number  £ for  n, /n  =0.1,  a/R„=0.05,  R~/R  =0.5,  and  several  values  of  q. 
bp  U 0 c 

From  Fig.  3,  we  note  that  the  number  of  unstable  modes  increases  as 
the  self-field  strength  (as  measured  by  q)  is  increased.  For  example. 


for  q=2,  the  mode  numbers  from  £-2  to  £=1  are  unstable.  On  the  other 


hand,  for  q=0.25,  only  the  £=4  mode  is  unstable.  Moreover,  the  maximum 
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growth  rate  is  increased  by  increasing  the  value  of  q. 

We  conclude  this  section  by  emphasizing  that  the  growth  rate 

exhibits  a sensitive  dependence  on  the  location  of  the  '■'-'•'.ducting 

wall.  This  is  illustrated  in  Fig.  4 where  the  normalized  growth  rate 

Imio/w  , and  the  normalized  Doppler-shifted  real  frequency  Refoi+Ho)  ,)/w  . 
C1  ci  ci 

are  plotted  versus  Rq/Rc  for  1= 3,  nb/np=0.2,  a/RQ=0.05,  and  several 
values  of  q.  Evidently,  the  growth  rates  as  well  as  the  Doppler  shifted 
real  frequencies  are  substantially  reduced  whenever  the  conducting  wall 
is  located  sufficiently  close  to  the  P-layer. 
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V.  CONCLUSIONS 

In  this  paper,  we  have  investigated  the  electrostatic  stability 
properties  of  a rotating  P-layer  immersed  in  a uniform  background 
plasma.  The  equilibrium  and  stability  analysis  (Secs.  II-IV) 
was  carried  out  within  the  framework  of  a hybrid  (Vlasov-fluid)  model 
in  which  the  layer  ions  are  described  by  the  Vlasov  equation,  and  the 
layer  electrons  and  background  plasma  electrons  and  ions  are  described 
as  a macroscopic,  cold  fluid.  Moreover,  electrostatic  stability 
properties  were  calculated  for  the  choice  of  equilibrium  ion  distribution 
function  in  Eq.  (10).  Although  a slow  rotational  equilibrium  (Pq>0) 
is  stable  [see  discussion  following  Eq.  (54)],  it  is  found  [Sec.  IV] 
that  a high-density  fast  rotational  equilibrium  (P^<0)  is  unstable  for 
a broad  range  of  physical  parameters.  One  of  the  most  important 
conclusions  of  this  study  is  that  the  background  plasma  has  a large 
influence  on  stability  behavior.  In  particular,  the  instability  growth 
rate  is  significantly  reduced  for  sufficiently  low  background  plasma 
density.  Finally,  we  conclude  that  the  growth  rate  increases  with 
increasing  self-magnetic  field  (as  measured  by  q) . Moreover,  the 
characteristic  growth  rate  for  £>2  is  a substantial  fraction  of  the  ion 
cyclotron  frequency.  The  fundamental  mode  (£=1)  , however,  is  found  to 
be  stable. 
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FIGURE  CAPTIONS 


Fig.  1 Equilibrium  conf iguration  and  coordinate  system. 

Fig.  2 Plots  of  normalized  growth  rate  Imw/u)^ versus  n^/np  obtained  from 

Eq.  (56)  for  1=2,  a/RQ=0.05,  Rq/Rc=0.5,  and  several  values  of  q. 

Fig.  3 Plots  of  normalized  growth  rate  Imw/a)^  versus  l obtained  from 

Eq.  (56)  for  n,/n  =0.1,  a/R  =0.05,  R^/R  =0.5,  and  several  values 
bp  U U c 

of  q. 


Fig.  4 Plots  of  (a)  normalized  growth  rate  Imw/w  , and  (b)  normalized 

ci 

real  frequency  Re (w+Jlw^) /ui  ^ versus  rq/rc  obtained  from  Eq.  (56) 
for  1=3,  n,/n  =0.2,  a/Rri=0.05,  and  several  values  of  q. 
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